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ANALYSIS OF STARVATION EFFECTS ON HYOROOYNAHIC LUBRICATION 
IN NONCONFORMING CONTACTS 
by David E. Brewe* and Bernard J. Haorock 
National Aeronautics and Space Administration 
LeMis Research Center 
Cleveland, Ohio 44135 
SUMMARY 

Numerical methods were used to determine the effects of lubricant 
starvation on the minimum film thickness under conditions of a hydrodynamic 
point contact. Starvation was effected by varying the fluid inlet level. 

Tne Reynolds boundary conditions were applied at the cavitation boundary and 
zero pressure was stipulated at the meniscus or inlet boundary. The 
analysis is considered valid for a range of speeds and loads for which 
thermal, piezoviscous, and deformation effects are negligible. It is 
applied to a wide range of geometries (i.e., from a ba11-on>plate 
configuration to a ball in a conforming groove). Seventy-four cases were 
used to numerically determine a minimum-film-thickness equation as a 
function of the ratio of dimensionless load to dimensionless speed for 
varying degrees of starvation. From this, a film reduction factor was 
determined as a function of the fluid inlet level. Further, a starved fully 
flooded ooundary was defined and an expression determining the onset of 
starvation was derived. As the degree of starvation was Increased, the 
minimum film thickness decreased gradually until the fluid inlet level 
became critical. Reducing the fluid inlet level still further led to a 
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sharp decrease In the* ininimum film thickness. An expression determining the 
critically starved fluid inlet level was derived. 

The changes In the Inlet pressure buildup due to changing the available 
lubricant supply are presented In the form of three>d1mens1onal Isometric 
plots and also In the form of contour plots. ' 

INTAOOUCTION 

The effect of starvation In a hydrodynamically lubricated conjunction 
can be studleo by systematically reducing the Inlet supply and ouservlng the 
resultant pressure distribution and film thickness. This starvation effect 
can have a significant role In the operation of machine elements. For 
example, roller-end wear due to roller skewing can be a critical problem for 
high-speed cylindrical roller bearings. It Is desirable that the hydro- 
dynamic film generated between the roller end and the guide flange provide 
stiffness and damping to limit the amplitude of the roller skewing motion. 
However, at high rotational speeds the roller end and the flange are often 
subjected to a depletion in the lubricant supply due to centrifugal 
effects. In such cases, tne minute amount of lubricant available at the 
roller-end-flange conjunction might well represent an example of steady- 
state starvation. Starvation effects in hydrodynamically lubricated 
contacts are Important also If one wishes to calculate the rolling and 
sliding resistance and/or traction encountered In ball and roller bearings 
Iref. 1). in another example, the effect of restricting the lubricant to a 
roller bearing is seen experimentally and theoretically to reduce the amount 
of cage and roller slip (ref. 2). The theoretical analysis was accomplished 
by Changing the location of the boundary where the pressure begins to build 
up and noting the effect on the hydrodynamic forces. Combining this with 
relative velocity expressions and equilibrium equations enabled the 
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determination of the amount of cage and roller slip. 

The location of the inlet and exit boundaries as well as the respective 
boundary conditions to be applied has been one of the most controversial 
issues concerning starvation of hydrodynamic contacts. The issue of the 
effect of the lubricant supply on the inlet boundary condition and its 
consequences to incipient pressure buildup began to materialize as a result 
of earlier studies applied to rigid cylinders (refs. 3 and 4). Lauder (ref. 
5) and Tipei (ref. 6) asserted an upstream limit of the fluid film where the 
pressure begins to rise as governed by the Reynolds equation. This limit 
according to Lauder is determined by applying >^everse flow boundary condi- 
tions (i.e.. u = au/ay = 0). Tipei locates the upstream limit as defined by 
tne line of centers of two bounded vortices that are observed for pure 
rolling. Both cases, nave been criticizeo because their analyses lead to one 
position of pressure buildup regardless of the oil supply (ref. 7). Oowson 
(ref. 8), Floberg (refs. 9, 10), and most oramatically. Wedeven, Evans, and 
Cameron (ref. 11) provide experimental evidence supporting the idea that the 
locatioh of incipieot pressure rise is determined by the oil supply. 

Further, wedeven, et al., using a Grubin type of EHO analysis, obtained very 
good correlation between experiment and the theory of starvation effects by 
choosing the start of the pressure buildup to occur at the meniscus 
boundary. Oteri (ref. 12), using stream function analysis for rolling rigid 
cylinders, showed that incipient pressure rise occurs at the meniscus 
boundary even in the presence of reverse-flow conditions. In view of this 
work, starvation effects in machine element applications can be predicted 
and relied on with a greater degree of confidence. 

One of the more important manifestations of lubricant starvation is the 
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reduction in film thickness. This topic has recleved a good deal of atten- 
tion In the literature (refs. 11, 13-21). With the exception of references 
13 and 19 to 21« these references are applicable only to elastohydrodynamlc 
situations. Most of the work concerned with rigid contacts has been devoted 
to Mne contact applications (refs. 13, 20, and 21). Oalmaz and Godet (ref. 
19) analyze the effect of the Inlet on the film reduction factor for a 
spnere against a plate. However, to the authors' knowledge, an effort that 
parallels that of Hamrock and Oowson (refs. 14, 15) for the EHO contact Is 
absent from the rigid contact th> ory. In those works, an ex|.>ress1on was 
determined that relates the tilni reduction to the inlet distance. 

The current study Is a resumption of a previous rigid-contact analysis 

(ref. 22) to extend validity for the minimum-film-thickness equation derived 

there over a wider range of film thicknesses as well as to Include the 

effects of starvation In this equation. The start of the pressure buildup 

as determined by the Reynolds equation Is assumed to occur at the Inlet 

meniscus. The location of the cavitation boundary was determined by 

applying the Reynolds boundary conditions as discussed in previous work 

(ref. 22). The study applies to a wide range of geometries (I.e., from a 

ball-on-plate configuration to a ball In a conforming groove). Seventy-four 

cases were used to numerically determine (1) an equation relating minimum 

film thickness with the fluid Inlet level as well as with the dimensionless 

loao-speed ratio and geometry, (2) an equation predicting the onset of 

starvation, and (3) an equation predicting the onset of a critically starved 

conjunction. The resulting equations are valid for dimensionless minimum 

-5 -3 

film thicknesses Hq ranging from 5.0x10 to l.OxlO . Further, 
contour Isobar plots and three-dimensional Isometric pressure plots are 
presented. 
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bYMBas 


least- squares coefficients 
difference, [(Hq - Hq)/Hq)] x 100, percent 
dimensionless film thickness, h/R^ 
dimensionless minimum (central) film thickness, 

dimensionless calculated minimum (central) film thickness 
dimensionless fluid Inlet level 

dimensionless fluid Inlet level (onset of starvation) 
film thickness, m 

minimum (central) film thickness, m 
reduced hydrodynamic lift, dimensionless 
direction normal to boundary 
dimensionless pressure, pR^/hqU 

pressure, N/m^ 

R R 

effective radius of curvature, w - 4. - ^ ■, m 

X y 

separation due to geometry of solids, m 

ratio of dimensionless load to dimensionless speed 

average surface velocity in x direction, (u^ ♦ Uy)/2, m/sec 

load capacity, N 

dimensionless coordinate, x/R^ 

coordinate along rolling direction, m 

dimensionless coordinate, y/R^ 

coordinate transverse to rolling direction, m 

radius ratio, RJR„ 

y ^ 

film reduction factor 


•’0 


fluid viscosity at standard temperature and pressure, N sec/m? 


0 

Suoscripts: 

A 
ti 

cr 
f 

x.y 

NUMERICAL PROCEDURE 

The hydrodynamic effects on the central film thickness between two 
rigid solids In lubricated rolling and/or sliding contact are analyzed under 
conditions of lubricant starvation. The effects of starvation are 
oetermlned oy systematically decreasing the fluid Inlet level. The Reynolds 
oounoary conditions are applied at the cavitation boundary, and zero 
pressure is stipulated at the meniscus or inlet boundary. The lubricant is 
assumed to be an Incompressible Newtonian fluid under laminar, isothermal, 
isoviscous, and steady-state conditions. The numerical approach follows 
tnai of a previous investigation (ref. 22). There, a fully flooded film 
profile was specified and a pressure distribution satisfying the Reynolds 
equation was determined for a given speed, viscosity, and geometry. The 
analysis treats the two rigid bodies as having parallel principal axes of 
inertia. This enables one to make a simplifying transformation to an 
equivalent system of a rigid solid near a plane separated by a lubricant 
film (fig. 1). 


Archard-Cowking side-leakage factor. 


m 


solid A 
solid B 
critical 

flooded conjunction 
coordinate direction 
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Relevant Equatiuns 

rne same dimensionless expressions are used here as In reference 22, 
that Is, 

X - x/R^, Y - y/R^. H . h/R^ U) 

also 

P “ pR^/n^u, d » 

Tne Reynolds equation 



Is the governing equation within the conjunction. 

We recognize that, when the Inlet supply levels are increased to values 
much greater than tne minimum film thickness, calculations as governed by 
tne Reynolds equation are inherently In error far from the center of 
contact. The reason Is that the Reynolds equation neglects curvature of the 
fluid film. Uowson (ref. 2i) has pointed out that the errors Involved In 
using this equation to determine the buildup of pressure in such regions are 
negligiole. Ihe predicteo pressures are themselves so very much smaller 
than tne effective load-carrying pressures in tne region of closest approach 
of tne solids. The dimensionless film thickness equation is given as 

H - ♦ 1 - o[l - - (Y/a)^j (i) 

where H is bounded above by the dimensionless fluid inlet level 

ana oelow by the dimensionless minimuiii film thickness (i.e., Hg £ H 
< e.g., fig. 2). The rluid inlet level is made dimensionless 

so that, it the irJet is completely filled, « 1 ana the conjunction Is 
said to oe fully flooded. 
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The Reynolds oounaary conditions are used* that Is* ? » • 0 at 

the cavitation boundary and P « 0 at the inlet boundary (H • A 

pressure distribution that satisfies these boundary conditions is then 
aetermined numerically by finite differencing. Taking advantage oi' the fact 
that the pressure distribution must be symmetric about the direction of 
motion* only half of tne domain was used in the calculations and the other 
half was aetermined from symmetry. A variable-mesh nooal structure was used 
to provide close spacing in and around the pressure peak. This helped to 
minimize the errors that can occur because of large gradients in the 
high-pressure region. As before* the integration of the pressure 
distrioution can be used in relating the hydrodynamic effects (i.e.* load* 
speed, and viscosity) to the minimum (central) film thickness for a given 
fluid inlet level. In general* for a given Hq and u tne load capacity 
w and/or the dimensionless load-speed ratio is determined as follows: 


W a OnUK 
0 X 


// 


P dXdY 


or 


where w/U « w/nyUfi^. 


w/u 


// 


P oXdY 


(4) 


fully Flooded Film Thickness 

The right side of equation (4) has been shown to be a function of the 
geometry and the minimum film thickness. For a fully flooded situation 
(H.^ s 1), equation (4) led to the following relationship in reference Zi 
(eq. (^0)): 

H,. « 128 aibUL/W)*^ for 10"^ < % < i^) 


where 


L - O.lil tan"^(a/2) + 1.683 


and 
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A Ulsaavantage of equation (b) Is that it cannot he applied to the full 
range of dimensionless film thicknesses normally encountered in a rigid 
— b 3 

contact (i.e.t 10 < < lU ). Consider* for example* a hall in 

rolling motion that is loaded against a flat plate* For a load-speed ratio 

M/U of 340* tne numerically determined value (i.e.* as determined hy finite 

-3 

difference analysis) of the dimensionless film thickness is 10 (table 
1). Equation (5) predicts the dimensionless ^ilm thickness to be 

1.22x10"'^, which is in error of the numerical value by 22 percent. 
Consequently, we wish to revise equation (5) so that it is valid for the 

thicker films as well. This revised equation should reduce to equation (6) 

in tne limit for tnin films. Furthermore* the revised dimensionless film 

thickness should be expressed in such a way as to easily include the effects 

of starvation. Tnis would then enable us to present one general expression 

to oe presented for the dimensionless film thickness that can be used for 

the full range of film thicknesses for a starved conjunction as well as for 

a fully flooded conjunction. After tne numerical analysis for each case was 

complete, the several curve fits or regression curves of on W/U 

were considered that would be consistent with these aoove requirements. The 

most suitable curve fit considered was in the form of a more general linear 

equation (5), that is. 
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A linear regression by the method of least squares (ref. 24) was 
performed by using the data (a > l.uu) in table 1 for the fully flooded 


cor.altlon to determine the constants In equation (6). The coefficient of 
determination r^ was found to be The value of r^ reflects 

tne fit of the data to the resulting equation: 1 being a perfect fit, and 0 
being the worst possible fit. The values of Cq and are given In 
table II. It can be verified tnat In going from eqiatlon (6) to equation 
(6), the coefficient of the load-speed ratio remains uncMnged and Is 
entirely a function of qeonietry, that Is, 

^ — 177 


When considering conditions in which the film Is thin. It can be shown that 
is insignificant to the dimensionless loau-speeo ratio term, whicn 
explains wny equation (b) worked so well for tne thin-filro regime. Using 
the coefficients C^, ana as oetermineu above, we car< determine the 
dimensionless film thickness from equation (6) as 



/ — ♦ 2.6bll 

\ bL yi23a ^ 


(S) 


This expression is valid fo** the full range of dimensionless film 
tnicknesses normally encountered in a fully flooded undefonned contact. 

Generalized Film Tnickness Formula (Applicable to 
biarved as well as Fully Floooed Conjuctions) 

Tne effect ot lubricant starvation on the hydrodynamic film thickness 
was ooserveo by varying tne fluid inlet level to the contact ana noting the 
effect on load capacity for five different film thicknesses for the 
ba I l-on-p late contact (i.e., a • I.UU). In addition to the fully flooded 
data, bb computer-generated data points were used to arrive at a family of 
equations having the form given in equation (6). An equation for each fluid 
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Inlet level was aetermlnea Dy pertorming a linear regression Oy the method 
ut least squares. Taole 11 lists for each fluid Inlet level the values of 
tne coefficients Cj and the coefficient of determination r^, anu 
1^..^ • the majuiiium percentage ut error U deftneu as 

HiqA 

u - ^ v X lao (9) 

% 

Tne taole Indicates that for all practical purposes* the constant 
is again strictly a function of geometry. Furthermore* the entire 
effect of starvation Is represented by the value of Cq in equation (6). 

An expression for the coefficient as a function of the fluid Inlet 

level would enable a determination of a generalized minimum-f 11m- 

thicKness formula that applies to starved as well as fully flooded 
conditions, a close examination of the variation of with in 

taole II reveals that Cy « ^or the severely starved 

situations. As aoproaches 1* Cy approaches a value very nearly 

equal to e (i.e.* the base of the natural system of logarithiRS, 2.7i6), 

H 1 n 

rms suggests e as a modulating factor. Further considerations for 
tne nearly flooded inlet levels snow that 

Co* 

Finally for the full range of values* varies with as follows: 



where x is a proportionality constant, equation (lu) can be improved by 
uoing an exponential curve fit using kite metnod of least squares on the 
data with H^„. This results in the 

r oi lowiitg equat lun: 
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( 11 ) 


l.lii/ e' 


U.tt604 H 


1/2 


The regression coefficient was 0.96U6« inoicating a reasonaolj^ good fit. By 
substituting (eq. (7)) and (eq. (11)) into equation (6) the 
minimum film thickness formula becomes 


w/U 


ftLUc'da) 


T77 



U.86U4 H 
e 
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( 12 ) 


in the interest of simplifying the equation still further, the curve 
fit constants were rounded off. The round off (i.e., 1.1U7 to 1.11, and 
0.B0U4 tu l.u) in this instance actually improved the agreement between 
and the actual in most cases. However, the round off results 
in a Curve fit that is no longer a minimum least squares deviation. Tl.'us, 
our generalized minimum film tnickness formula in terms of geometry (i.e., 
radius ratio a), loau-speeu ratio w/d, and fluid inlet level can be 

written as 


% “ 


ri/U 


bL( 128(1 ) 
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The measure of agreement between the calculated and input values of 
IS represented by tne value of u (eq. (^)) anu presented in table 111* 
Table lii snows tnat, when tne fluid inlet level is of tne same order of 
magnitude as the minimum film tn' ness, the error that results from using 
equation (U) oecomes larger. However, equation (li) can confidently oe 
used ror tne full range of minimum film thickness if the fluid inlet level 
IS suen tnat u.uoa < H^^ < i.uou. for very thin films (i.e., Hq < 

1(T^) equation (li) can be useful throughout the full range of fluid inlet 
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levels that were investigated (i.e., U.uOl < < l.OOU). Note that the 

film thickness formula is intended to be used only for a range of speeds and 
loads in which piezoviscous and deformation effects are negligible. It was 
determined that these effects can be significant for minimum film 
thicKnesses less than 5.UxiO~^. 

Note also from table 111 that excellent agreement is obtained for the 

near- line-contact applications (i.e., a » 36.54) even though these data were 

not used in the determination of the above equation (13). riost of the 

predictions by equation (13) are within 3 percent of the numerically 

determined values and do not exceed 6 percent for any case. 

$ 

Equation (13) is equally valid for the fully flooded conjunction as 
well as for the starved conjunction. That is, if we set s 1.00, 

then Hq - HQ^f and equation (13) reduces to 

V< • •■*) '' 

In comparing this with equation (8), one discrepancy appears. The value of 
the constant from this equation is slightly different than that determined 
for equation (8). This comes about, for the most part, from the 
simplification of the generalized expression (eq. (12)). This difference 
has not introduced any appreciable error. Table III shows that the error 
for the full range of minimum film thicknesses investigated (for * 1) 
is less than 1 percent for a « 1. for the near- line-contact geometry 
(i.e., a = 36.54) the error does not exceed 3.36 percent. 

Keduction in Minimum film Thickness 

It is now possible to determine the reduction in minimum film thickness 
from the fully flooded value if the fluid inlet level is known. This can be 
done by inserting equation (14) into equation (13), which gives 
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Uividing both sides of the equation by gives 



wnere 3 is the reduction in minimum fi'i.H thickness due to starvation. 

KESULTb ANU UISCUSSION 

effect of starvation on Pressure Distruoution 

Tne discussion of lubricant starvation can be facilitated by focusing 

on one of the simplest geometric arrangements (i.e., a ball rolling and/or 

sliding against a flat plate) as shown in figure 2. Tne figure compares the 

pressure distribution determined numerically for the ful ly-f loodeo inlet 

with tne most severely starved inlet = o.OOl). The comparison is 

made for a constant minimum film thickness (i.e., Hq = 1.0x10"^). Note 

that tne pressure peak built up in the starved inlet is only slightly 

smaller than that of the fully flooded inlet. However, the area of pressure 

0 U 1 Id-up is consioeraoly smaller, and so the starved inlet is unable to 

support as much load for a given film thickness as the fully-flooded inlet. 

figure i provides the same sort of comparison but for a thicker minimum 

film (i.e., Hq = 1.0x10”'*). The significant difference between the two 

figures is tnat the starved inlet for the tnicxer film has a more pronounced 

effect on tne pressure peak. The fluio inlet level (H, = U.UU2) for the 

in 

thicker film represents a relatively more highly starved inlet since 
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is of the order of in this case. The other feature to oe noticed in 
cotnparing figures 2 . and 3 is that the pressure distriDution is more evenly 
spread out tor the thicker tiitii. Thus, cnanges in the meniscus (or 
integration domainj are going to nave a more noticeable effect on the load- 
carrying capacity. Note also that because of the boundary conditions the 
integration oomain takes on a “kidney-shaped" appearance. This is more 
clearly shown in the isobar 1c contour plot shown in figure 4. 

Minimum Film Thickness Equation 

Thus far we have compared the pressure build up ih a severely starved 
inlet with that in a fully flooded inlet for two minimum film thicknesses. 
Our investigation, however, included several fluid inlet levels for a 
variety of minimum film thiCKnesses. The results are summarized in table 
Hi. 

A generalized minimum film thickness formula (eq. (13)) was derived 
from tne results of table III. Figure 5 indicates how well the equation 
represents tne computer generated data in tne table. It was not possible to 
display all these results in figure 5. However, the figure is 
representative of the overall results. Tne equation fits the data quite 
well except when the fluid inlet level is of tne same order of magnitude as 
the minimum film thickness. Based on the discussion concerning peak 
pressure, it would seem that the formula holds well for those cases in which 
tne degree of starvation is such that peak pressure is not signif icantly 
reduced. 

df course, it would oe most desirable to compare the data in table III 
with experimental data. To tne authors' knowledge, the only available 
experimental data were obtained by ualmaz and Godet (ref. 2t>). To compare 
equation (13) with experiment, the data from reference 2b were replotted in 
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figure 6. The experimental aata were taken under lightly loaded (rigid 
contact), isoviscous conditions for pure sliding of a oall on a plate. The 
fluio inlet level in these experiments was reported to oe 1 millimeter. The 
oall oiameter was JU millimeters; consequently the dimensionless fluid inlet 
level was U.U67. The experimental data were presented as a plot of 

the dimensionless parameters H^/WG versus The 

materials parameter G in the plot was included so as to accomodate the 
elastohydrooynamic range in a more general way. 

Here, we wisn to compare our hydrodynamic starvation theory only with 
tne hyorodynamic results of Daimaz and Godet. To do this, the ordinate and 
doscissa were properly scalea (assuming a reasonable value of 
Wb = 4,oxlu“^) so that the minimum film thickness coulo be plotteo against 
tne dimensionless load-speed ratio as snown in figure 6. The solid 
line in figure 6 is a plot of equation (ii) for o = 1 ana » 0.067. 

The aasheu line represents a previous theory (ref. dZ) in which the reduced 
inlet level due to starvation was not consiuerea. The present theory (eq. 
(id)) snows better agreement witn the experimental results than our previous 
theory (ref. Z ) for lower values of W/U. This is because under conditions 
of low loau anu/or high speed, the pressures become more significant away 
from the center of the contact. Consequently, neglecting tne size of tne 
inlet introauces an increasing amount of error as M/U is decreased. 

Although consideration of the size of the inlet domain improves the 
agreement of tne theory with experiment for the thicxer films, still further 
improvement is possible. It is believea that the effects of reverse flow in 
tne inlet must oe incluueu in tne theory to obtain better agreement. If 
reverse flow is consiaered, not all the available lubricant determined by 
tne fluia inlet level will pass tnrougn the contact. The hydrodynamic 
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contact would essentially see tnis as a reduction in supply from what 
actually is there. In other words* the inlet is more severely starved than 
we nave taKen into account. From figure 6* we see that increasing the 
severity ot starvation has more effect on the load capacity for thicker 
films (i.e., than it does for thinner films (i.e.* « 

Consequently* reverse-flow considerations should improve agreement 
oetween theory ano experiment for the thicker films while still maintaining 
good agreement in the thin-film range. 

LuDricant Film Thickness Keduction Factor ano Onset of Starvation 

Of practical importance to lubricant starvation is the reduction in 
minimum film tnickness s from the fully flooded value. Equation (16) is a 
derived expression for a in terms of the fluid inlet level and the fully 
flooded film thickness (also given oy eq. (14)). Figure 7 is a olot of a 
as a function of tne fluio inlet level for several values of 

Ho,f* It is of interest to determine a fully flooded - starved Doundary 

(i.e.* tnat fluid inlet level after which any further decrease causes a 
significant reduction in tne film thickness), riamrock and Uowson (ref. 14) 
determined tnis boundary for elastohyarodynamic (tHU) applications upon 
satisfying the following condition: 

I in in 

fne value of U.U3 was used in equation (17) since it was ascertained that 
the data in taoie I were accurate to only percent. 

Thus* for a given value of one can solve for a value of 

* 

tnat satisfies equation (17). A suitable relationship between 

and can be obtained oy generating a table of values 

(e.g.* table IV) ano fitting a power curve by the method of least squares. 
This gives 
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Thus, we nave an equation tnat determines the onset of starvation. That Is, 
* 

for H. > a fully flooueo condition exists, wnereas for 
in — in 

^1n ^tn ® starved condition exists. 

Critically Starved Inlet 

In certain bearing applications, the power loss resulting from churning 
of the oil may be higher than the power loss resulting from friction of the 
bearing alone (ref. ^6). These power losses can be minimized by reducing 
the lubricant supply until a loss In film thickness causes the friction 
losses to Increase. According to the results shown In figure 7, the fluid 
Inlet level can be decreased substantial ly without adversely affecting the 
minlmuiii film thickness. Consequently, It might be advantageous to operate 
cne bearing with a lubricant supply just sufficient to preclude any drastic 
reductions in minimum film thickness, as seen In the figure. Such a 
critical fluid Inlet level might well oe defined to occur at the knee of the 
curve, that Is, 

= 1 m 

sH ' 

in in,cr 

Solutions to this expression for several values of minimum film thickness 
are listed in table IV. A power curve fit oy the method of least squares 
gives 

"in.cr - 

The regression coefficient for this fit was determined to be 0.9^9^, 
inoicating an extremely good fit. 
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CUNCLUUlNti REMARKS 

Numerical methods Mere used to determine the effects of luoricant 
starvation on the minimum film thickness under conditions of hydrodynamic 
point contact. Starvation was effected by varying the fluid Inlet level. 
The Reynolds boundary conditions were applied at the cavitation boundary, 
and zero pressure was stipulated at the meniscus or Inlet ooundary. The 
analysis Is considered valid for a range of speeds and loads for which 
thermal, plezoviscous, and deformation effects are negligible. It can be 
applied to a wide range of geometries (l.e., from a ball-on-plate 
configuration to a ball In a conforming groove). Seventy-four cases were 
used to numerically determine 

(1) A generalized expression for the minimum film thickness as a 

function of dimensionless load-speed ratio, geometry, and fluid Inlet level 

(eq. (13)). The expression should be applied for film thicknesses In the 

range i.Oxlu"^ i i 1.0x10“^ and for fluid Inlet levels of U.004 < 

-h -4 

< 1.00. For 5x10 1 rij < Id the equation can be applied for 

a fluid Inlet level of 0.001 < < l.oo. 

( 2 ) A film thickness reduction factor (eq. (16)) expressed as a 

function of tne degree of starvation (or fluid Inlet level) for a given 

fully flooded film tnickness value. 

(3) An equation (eq. (Id)) that determines the onset of starvation. 

(4) An equation (eq. (^0)) that determines a critically starved 

contact. Contour Isooar plots ana three-dimensional Isometric plots are 




also presentee. 
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IAtM.£ 1. 

RATIO FOR FIlTE VALOFS OF MININUM FIlH TH|CANt&& 


[Avcrtge surfAcr veloctty tn a-direction, u, 10 ctWicc; fluid vitcotlty. 
0.41U1U~^ N %K/ci»i effective reoiui of curvature, ft,, l.ilU^ cn.J 


Ulmens ion less 
f luid Inlet 

level. 

Olmenslonless minimum film tnicaness, Hj 

UlO*^ 

b»10*!* 

1> 

10*^ 

b»10'* 

lalO-'* 

*l|n 



Ha^ us 

ratio, a 




1.0 

1.0 

l.O 

Oo.bO 

1.0 1 1.0 



Dimensionless load-speed ratio 

N/0 


1.000 

./HO 

.bOU 

.diO 

AM 

.0/0 

.00b 

.0/0 

.010 

.004 

.00/ 

.001 

J/00.14 
J/0b.b4 
i/OJ./d 
ObRb.a/ 
Jotib.41 
iODb.40 
i04ti .44 
OolO.Ul 
ibb/.bb 
044/. eJ 
0/U4.44 
Olod.bd 

1640. /4 
1646.// 
1040.41 
loOb.Otf 
lo/O.Ob 
loOo.lo 
Ib/O.ol 
IbOo. JO 
14/4.64 
10//. 16 
l//l.b/ 
1041. bO 

llbO.64 
llb0.O4 
llbO.oO 
114/. 60 
110/. 6/ 
1106.40 
10//. 46 
1044.01 
46/. 60 
66/. b/ 
/Ob.ol 
bo/./b 

1/400.4/ 
1/4/6.06 
1/400.// 
1/0/0.66 
1///4.01 
1/04/. 0/ 
llbb/.l6 
11164. /O 
10b00.40 
4/0/ .64 
/604.46 
61//.10 

44b. b4 
44b. 00 
44/. 6/ 

464.66 
4/4. /O 

460.66 
4/1.74 
067. /O 
00/. 44 
/0/.77 
1//.47 

07.06 

004.6/ 

004.07 

006.64 

0/6.74 

018.46 
Z46.66 
/6b. 76 
/04.64 

164.46 
64.74 
/b.4b 
44.60 


*fnis value of u/U detemined Dy usinq • F.iiilO'*, 
since Hq • UlO*-^ results in no load capacity. 
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lal Fully flooded condition; dimensionless fluid inlet level Hjp ■ 1. OQi 



(b) Starved condition: dimensionless fluid inlet level H|p*0i00l. 

Figure Z. * Three'dimensional representation of pressure distribution, comparing starved with fully 
flooded conjunction for dimensionless minimum film thickness Hg • 1x10 . 




!3) Fully flooded condition: dimensionless fluid inlet level H- ■ l.OO. 



starved condition.- dimensionless fluid inlet level H,-, - o. OOZ 



(a) Fully flooded condition: 
dimensionless fluid inlet level 
Hjn • 1.00; dimensionless 
maxim^ pressure • 
1.20x10°; dimensionless load- 
speed ratio W/U ■ 1153.6. 


(b) Starved condition: (c) Starved condition: 

Hjn ■ a 004; P^ax * ^^in * ^ ^max * 

W/U-86^6. W/U- 567.8. 


Figure 4. - Isobaric contour plots for three fluid inlet levels for dimensionless minimum film thickness 
Hq ■ 1x10'^. 
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Eq. (13) 



Dimensionless load-speed ratio, W/U 


Figure 5. - Comparison of dimensionless minimum film thickness equation leg. (13)1 
with computer-generated data as a function of dimensionless load-speed ratio for 
several values of dimensionless fluid inlet level. 
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Figure 6. - Theoretical and experimental results. 
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Figure 7. - Minimum film thickness reduction factor as a 
function of fluid inlet level for several values of dimen- 
sionless minimum film thickness for flooded conjunction 



